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Abstract. This paper is devoted to guaranteed 
estimation 1 of linear functions, defined on the so- 
lutions domain of the linear descriptor difference 
equations (LDDE) system, where right-hand part 
and initial condition are arbitrary elements of the 
given set. Minimax estimations are build on the 
basis of system's state observation with unknown 
deterministic noise. Minimax filtration task is 
studied for LDDE system with special structure. 
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Introduction. 

Suppose that vector {x k }^ +1 satisfies linear de- 
scriptor difference equation 



Fk+i x k+i ~ C k x k — Bkfk, k — 0, N 
with initial condition 

F x = 



(2) 



where F k , C k - m x n-matrixes, B k is mxp-matrix. 
We'll be interested in building minimax approxi- 
mation of the linear function 2 



N+l 



{{ x k}) = y~l(4,fffc)n,4 e 



fc=0 



assuming that 

HI state Xk observations are given in the form 

of 

y k = H k x k + g k , k = 0,N, (3) 



where {g k }o is some deterministic noise,H k is q x 
n matrix and 

H2 {/fc}^, {gk}o ^re some arbitrary elements of 
the ellipsoid 

G ={({/*}*, {flfc}S r ):(g- 1 /-i,/-i) P + 

A (4) 

/ XQkfki fk)p + (Rk9k,9k)q ^ 1}, 
k=0 

where Q k , R k are some symmetric positive-defined 
matrixes with appropriate dimensions. 

Assume that motion of some object 3 is de- 
scribed by LDDE (0) with initial point that satis- 
fies (J2J) while system disturbance (f k ) along with 
right part in (J2J) and noise (g k ) in the object's 
state observation model (jlj) are supposed 4 to be 
uncertain. Than mentioned above problem can be 
treated as guaranteed estimation of the object's 
transfer (from the set of possible initial states de- 
scribed by (J2J), caused by uncertain disturbances, 
on the basis of noisy observations. Among an- 
other applications of the guaranteed estimation 
task studied in this paper there is a image mod- 
elling [1^ and constrained robots movement |2j. 

Let us introduce a notion of minimax a- 
posteriori set. At first we'll define set Af clS db col- 
lection of all pairs ({£/c}o + \ {fk}^) satisfying ((T}- 
(j2J). Than, let us set 

N+l 

JMhY-v {^-}o +1 ) = £(Q*-iA-i, fk-i) P + 



k=0 



(Rk(Uk - H k x k ),y k - H k x k ) 



(5) 



where y N+1 = 0, H N+1 = 0. If ({f k }% G G 

satisfies (JI}-© for some {^fc}o +1 , than 

({fk}^, fen eAr,j v ({f k r_ v {x k r ^) < i 

(*) 

and vice-versa if ({fk}^, { x k}o +1 ) satisfies (*) 
than 



({fk}-i,{9k = Vk 



H k x k } 



1 So-called minimax estimation 
2 (-, •)„ denotes inner product in 



3 Lot's of examples we can find in robototechnics 0E| 
4 For instance, ({fk}- lt {gk}o) could be measured only 
with some errors; ({/fc}^, {gk}™ ) they are random but we 
do not have exact information about corresponding corre- 
lation functions. 
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It means that we can describe a set of all {x k }^ +1 
causing to appearance of given {y k }^ in (jSJ) while 
({/fe} N n idk}o) run through some subset 5 of G. 



Definition 1. The collection 



({f k }^,{yk-H k x k }^)eG} 
is called a-posteriori set. 



(6) 



It's obvious that real solution {x k }™ +1 of CEJ)-© 
being observed in Q for some ({/fc}^, {gk}o ) £ G 
belongs to G y . Hence it's naturally to look for the 
i({xk}) estimation only among the numbers from 

L={t({x k })\{x k }^£G y } 

Because of uncertain ({/fc}^, {gk}o ) we'll use 
minimax strategy for finding optimal estimation 

£({x k }) from within L: for each {x k }^ +1 E G y 
we need to calculate greatest distance between 
£({x k }) and L - so called guaranteed estimation 

error a 

£({x k }), where {x k }™ +1 has a minimal a({x k }^ +1 ). 

Definition 2. Linear function £({x k }) is called 
minimax a-posteriori estimation if 

inf sup \£({x k })-£({x k })\ = 

{x k }£G y { Xk }<zG y 

sup \£({x k }) - £{^ k })\ 

{x k }GGy 

The non-negative number 

a= sup \£({x k }) - £({x k })\ 

{x k }£Gy 

is called minimax a-posteriori error. 

In next section we shall study criteria of mini- 
max a-posteriori estimation existence along with 
minimax error finiteness. It'll also be discussed a 
few ways of minimax estimation calculation. 



5 This subset consists of all pairs ({fk}-i, {.Sfcfo ) e G 
satisfying f° r some {xk}„ +1 - 



Minimax a-posteriori estima- 
tion. 

Theorema 1. If £({x k }) = X]fc=o(4, x k ) n and 

{4} +1 e L = {4 = F' k z k - C' k z k+l + H' k u k , 
F^ +1 z N+1 = £ N+1 , z k e M m , u k e W} 

than 

N+l 

=X)(4,£*)n, (7) 



k=0 
N 



N+l 



o- = (l - ^(yk,Rk(yk - H k x k ) q ) 2 Q2(£ k ,Pk))' 



k=0 



k=0 



(8) 



where {x k }^ +1 is a solution of 

F k z k ~ C' k z k+1 = H' k R k (y k - H k x k ), k 
F k +ix k+ i — C k x k = B k Q k B k z k +i, 
F x = B-\QZiB'_ 1 zq 1 Fn +1 zn+i = 



0,N 



(9) 



and {p k }o +1 is a solution of 

F k+1 p k+1 = C k p k + B k Ql l B' k d k+ i, k = 0, N, 

F k d k = C k d k+ \ + £k — H' k R k H k p k , 

-Fn+i^n+i — ^N+i) ^oPo — -B-iQ-i-B-i^o 

Next theorem gives a recurrence algorithm for 
minimax a-posteriori estimation calculation in 
case of special structure of matrixes Fk,Hk. We 
also suppose here that number of measurements 
is equal to N + 1 hence it's not necessary to set 
y N+1 = 0,# N+1 = in ©. 

Theorema 2. //rank ^ = n and Bk = E than 
for any £ G W 1 minimax a-posteriori estimation 

(£, xn) of inner product (£, x^s) n can be represented 

as 

(I -r x ) = (£,x NjN ) n 

where 

&k\k = Pk\kF' k (Q k \ + Ck-lPk-l\k-lC' k _i) X 

x Ck-\Xk-i\k-i + Pk\kH k Rkyk, 

Pk\k = [E k {Q k *x + Cfc_iPfc-l|fc-lCfe-i) 1 Fk + 



H' k R k Hk) , Po\o = (FqQF + H' R Hq 
Xq\o = Po\oHQR y i 



.-I 



(10) 
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Example. 



Proofs. 



We'll illustrate theorem |2] in case of estimating 
inner product (a, Xn+i) for linear stationary de- 
scriptor difference equation. Let us set m = 2, n = 

3,/ = 4,N = 200, /_! = F k 
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and let's choose {fk}^ and {gk}o from the unit 
sphere. Simulated state {xk}g +1 and it's minimax 
a-posteriori estimation {xk}^ +1 are shown on the 
figure H 
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Figure 1: {xk}„ +1 (blue, dashed) and {xk}g 
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Theorem 's proof. The following lemma gives criteria for guaran- 
teed estimation error finiteness. 

Lemma 1. 

sup \£({x k }) - £({x k }™ +1 )\ < +00 & {4}* +1 6 L 
{x k }eG v 

It's easy to see that sup^^^, \d — c| = D + [c| for any real D, c. 
This implies to 

sup \£({x k }) - £({x k }™ +1 )\ = 

i[ S ({4} N+1 |G s )+ S (-{4}? +1 |G«)]+ (11) 
K({2fe}? +1 ) - ^[s({4}o +1 |G„) - s(-{4}? +1 |G H )]| 
for {4}o +1 6 L, so 

<(fo}) = i[s({^}? +1 |G w )-s(-{4}? +1 |G !/ )], 
&=^s({4}? +1 |G J/ )+ S (-{4}? +1 |G a )], 
where s( -|G H ) is a support function of G y . Let's find s( -\G y ). 

Lemma 2. Vector ({fk}, { x k}o +1 )> w here f k = B' k z k+ i, k = — 1, N 
and z k ,x k are some solutions of J<JJ, is a minimum point of the J y 
on N 



(12) 



: def 



N 



J"= inf Jy = V" (RkVk,Vk ~ H k Xk) q 



and 



Jy({fk - fk^Axk -Xk}) = Jb({/*}? 1 ,WJ +1 ) + J, (13) 
/or W 5 ({/fe}"i, {^fclo +1 ) 6 AT 

If we set P({/ fe } N 15 {x fc }o +1 ) = {^fclo +1 and 

G d ^{({/ fe }^, {x h }? +1 )|J,({A}!! I1 {x fc }? +1 ) ^ 1} 

than Gy = P(G PljV) as it follows from G H and G definitions. For- 
mula 1131 implies to 

GnA/" = ({fk}, {xk}o +1 ) + Go nyV, 

where 

Go : ^({/fc}-!, W>0 +1 ) < 1 " J} 

It's easy to see that Go = —Go- Hence for {4}o +1 G £ 
N+l 

s({4}o +1 |G y ) = £ (£ k ,XkU +s({4}o +1 l^(G n AO) 

fc=0 

Last formula with a regard to 1121 implies to 171 , 

To prove J^J we need to find support function of the set Go n 
J\f . This task |7J p. 164, c. 16.4.1] is equivalent to minimisation of Go- 
support function over affine set f ({4}) — N ± ■ It's easy to see that 

Argmfs F , ({< , fc}) _^_L C doms(-|Go) nP'({4}) - A^ 

Taking into account |jj p. 136, T. 13. 5] we can show that 

doms(.|Go) = { ^Itt - k = O^ + T} W 

so regarding to structure of L it's not difficult understand that min- 
imax error can be represented as JSJ - □ 



3 Jo = J y for {y k }» = 0. 
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Collorary 1. Suppose that all assumptions of the theoremTHare ful- References 
filled. To find a representation of the minimax a-posteriori estima- 
tion £({x k }) than it is sufficient to find a solution {{f k }, {xk}o +± ) 
of optimisation task 
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Proof. As it was mentioned above we can represent minimax a- 
post erior i estimation in terms of a-posteriori set's support func- 
tion 1121 . From the other hand each element of G y could be treated 
as sum of {£fc}o +1 and some vector from balanced set M n Go- □ 



Theorem proof. Theorem conditions implies L is equal to the 
whole Euclidean space (R n ) N + 1 , so (0, ...,£) £ L for any t e 1" and 
thus we can use corollary^from where and according to theorem 
conditions we need to find a solution of 



N+l 

J N +i({a:fc}o +1 ) = W F k x k ~ C k -lx k -l\\Q k _ 1 + | life - H k x k \\ 2 Rk , 

k-Q 

C_i = 0,x_! = 

(14) 

It is shown in |5] that we can obtain solution of 1141 using recurrence 
process Hot . □ 
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Lemma's 1 proof. Let {£ k }g +1 G L. It's easy to see that 



sup \£{{x k }) - <({i fc }? +1 )| < +co 
{x k }eG v 

SUp |^({Zfc})| < +oo. 



On the other hand if we set B—i = B,f-± = / than after some 
simple calculations we obtain 



^({zfc}) = ^2( F k z k ~ C' k z k+1 ,x k ) n + (H k x k ,u k ) q + 



(z N+ i,F N+ ia: N+ i) m = (y k - g k ,u k ) q 



N+l 



£( B *-iW*-i)p <+°°. V( It ); +1 eG s 



because of H2. 

Now we'll assume that {^fe}o +1 ^ L. It implies that 

£({xk}) + for some {x fc }g +1 : ({F k x k ,C k . 1 x k _ 1 , H k x k }) = 0, 

s°sup {xk}eGy \e({x k })\=+co. □ 



Lemma's 2 proof. Taking into account special structure of J y and 
vector's projection theorem in Hilbert space it's easy to prove that 

AT d = Arginfjv J y ^ and 



k = 



ifk-i, f k -i)p + (H' k R k {y k - H k x k ),x k ) n = 0, 



V({/ fc } N 1 ,{^}? +1 ) e N if «A},{i fc } N+1 ) e N. it implies 

(Qf,{Qkfk},{H' k R k {y k - H k x k )}) G Af x . On the other hand 
(/, {f k }, {x k }g +1 ) £ Af. Now with a regard to structure of J\f we 
can show that exists {z k } which satisfies JgJ and f k = B' k z k +i,k = 

^I7n. □ 
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